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The influence of chaos on properties of dilute instanton gas in quantum mechanics 
is studied. We demonstrate on the example of one-dimensional periodic potential that 
small perturbation leading to chaos squeezes instanton gas and increases the rate of 
instanton tunnelling. 
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1 Introduction 



Semiclassical properties of classically chaotic dynamical systems (quantum chaos [|Tj] ) is a 
reach rapidly developing field of research. One of the attractive phenomenon obtained in 
this direction is chaos assisted tunnelling (H and references therein). 

In this work we analytically investigate influence of the perturbation leading to chaos on 
the rate of instanton || transitions on the example of one-dimensional periodic potential. 
Instantons represent alternative way to describe quantum tunnelling and play an important 
role not only in quantum mechanics (recent example |4]]), but mainly in the modern theories 
of particle physics, e.g. electroweak theory, QCD, SUSY, etc. (for review see ||). The 
search of instantons, initiated in DESY, is one of targets of modern QCD ||. The main 
problem in this direction is an extremely small probability of instanton-induced events. 
It is of importance to study if chaos can increase probability of instanton tunnelling. We 
demonstrate on the particular example of simple quantum mechanical system that it is really 
possible. 

Hamiltonian of the system under consideration is taken in the form 

H = -p 2 + Uq cos x — ex S(t — nT), (1) 
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T is the real time period of perturbation. The systems with spatially periodic potential are 
well-studied in solid-state physics M and instanton physics M . Perturbation used in (|T]) was 
widely exploited in the systems exhibiting quantum chaos ||. 

For applying instanton technique we have to consider classical solutions of Hamilton 
equations in imaginary (Euclidian) time. Hamiltonian ([I]) has the same form (translated on 
7r) in Euclidian time as in real one. 

In Euclidian time classical Hamiltonian of the system looks as follows H = Hq + V, 

H = ip 2 - luI cos x, (2) 

+00 

V = ex S(r-nT). (3) 

n=— 00 

Here Hq is non-perturbed Hamiltonian of the system and V is the Euclidian potential of the 
perturbation. In these expressions variables without tilde denote corresponding quantities 
in imaginary time. 



2 Instanton gas in non-perturbed system 

Let us consider at first non-perturbed Hamiltonian ([I]), i.e. put e = 0. The system is 
characterized by degenerated vacuum structure on the classical level: 

x v n ac (t) = \n >= irn, n = ±l,±3,... (4) 

As is well-known any tunnelling transition in this system can be represented in terms of path 
integral in imaginary time: 

< m\e~ H ° r \n >= J[Dx]e~ s[x \ (5) 

where T is a time of the transition between wells, S[x] denotes Euclidean action. The main 
contribution to (|5|) is given by instantons. One-instanton configurations are classical solutions 
of Euclidean equations of motion and describe tunnelling between neighboring vacua. They 
can be easy found as well as one-instanton action (see, for example ||) 

x mst (r, t ) = ± (arctg [ e ^~^] - vr) , S[x mst (r, r )] = S inst = 8lu , (6) 

where arbitrary parameter To is a center of the instanton; sings '+' and '— ' correspond to 
instantons and anti-instantons. 

Multi-instantons are not exact classical solutions, but they give leading contribution to 
the amplitude of tunnelling fl5|) between distant wells 

n 

WiVi =£V n * t (T,r < ), 5[x (n) (r)] =nS mst . (7) 



X y 'IT 



t=l 



Here we suppose that time intervals between centers of single instantons 7* are not too close 
to each other (dilute instanton gas approximation) 0. 
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Let us focus our attention on the amplitude A q (T) of q one-instanton transitions dur- 
ing time interval T. The difference between instantons and anti-instantons is not essential 
in these calculations. Using standard instanton technique M it is easy to obtain Poisson 
distribution on q for the amplitude in Gauss approximation: 

A q (T) = N^e-^ (T^/S^) 9 , (8) 

where factor iV provides correct normalization, A (T) means absence of instantons. The 
average number of instantons for the time T is 

< q >= e- smst Y^/S^. (9) 

Thus in imaginary time we have a gas consisted of instantons (|6]) with average time interval 
between them i] and average density p : 

-p, ginst 

Vo = ~ L - = -^==, Po = Vo 1 = e~ sin3t VS^ = 2v^e- 8 -. (10) 
< q > VS 1 

Here, as usual, we neglected by the instanton size. 

Let us suppose large value of one-instanton action (§). It corresponds to high energy 
barriers. In this case we obtain strongly rarefied instanton gas. 

It should be noted that result fllTf) was obtained by using of both exact (|6|) and approxi- 
mate (0) classical solutions. Taking into account only exact solutions leads to zero instanton 
density p (or infinite interval rj ). Finite density can appear in non-perturbed system ([[]) 
only due to approximate solutions contribution. 



3 Squeezing of instanton gas due to small perturbation 

Now we consider the case with e ^ and estimate an average interval between instanton 
transitions (inverse density of instantons in instanton gas) for the perturbed system. For 
this purpose we represent the perturbation (Q) in the form 




V = -x\2 2^ cos(mz/r) + l . (11) 



Here v = 2tc/T. 

Perturbation (|TT|) destroys separatrix of non-perturbed system (||) and on its place 
stochastic layer appears PH| . We primarily estimate the width of stochastic layer. We 



use method reviewed in We apply it to the system (fj-|[) which was not considered in 
that works, although non-perturbed Hamiltonian (fj) and perturbation (|) were used inde- 
pendently from each other. 

Exact equation of motion for action variable is 



' lI -H = [ 2 £ cos (mur) + 1 ) , (12) 

m=0 / 



dHn ujT 



Here / denotes action variable and u = dHo/dl is a frequency of nonlinear oscillations (see 

H)- 
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Consider behavior of the system near separatrix (in imaginary time). Dependence of 
velocity x on time has the form of rare soliton-like impulses. Each impulse corresponds to 
rapid transition between two neighbor peaks of potential. Long time interval between two 
impulses is the time to get over a peak. It is seen from equation ( |12|) that action variable 
changes mainly during the impulse of velocity. Introduce phase of external force tp defined 
as up = v. 

To study chaotic behavior of the system we transform differential equation (12) to discrete 
mapping 



--I + 
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C{I,<p) 



(13) 



where (J, <p) denote values of action variable and phase just after the impulse of velocity, 
(J, uo) are the same quantities after previous impulse and 



r ( + °° 
C(I,<p) = - drxhY, 



cos (m<p(t)) + 1 



(14) 



Here we integrate over time interval of velocity's impulse At. From (|13|) we obtain (for small 

itv ev 2 duj 

+ (15) 

Therefore parameter of local instability (defined in ji~l]l ) is 



K 



5<p 
5<p 



1 



neis 



duo 



dl 



Co, Cq = 



dC 



dip 



(16) 



In the vicinity of separatrix Cq can be calculated 

sinh {ttv/2uq) 



C n ^ 47T 



cosh (7tu/2oj ) 



(17) 



Using (^) and estimation for \du/dl\ at \H — H s \ <C H s we can represent the parameter of 
local instability flT6"| ) in the form 



K 



2eu 2 sinh {ttu/2ujq 



lj cosh {nu /2ujq)\H — H s 



18) 



Here H s = uj 2 is the energy of non-perturbed system on the separatrix. 

Condition K > 1 means that dynamics of the system is locally unstable. Local instability 
leads to mixing and chaos | H| . Thus condition K — 1 gives the estimation for the width of 
the stochastic layer as follows 



\Hh — H K 



Aeu 2 



(19) 



under the assumption that v > ujq. Here is an estimated value of energy on boundary of 
the stochastic layer. A set of trajectories appearing in the stochastic layer can be considered 
as new multi-instanton-like exact solutions of Euclidian equations of motion for perturbed 
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system. These solutions demonstrate finite intervals between instanton transitions. The 
reason is that for majority of trajectories in stochastic layer finite time is needed to pass 
from one potential well to another. Thus density of instanton gas strongly increases in 
comparison with density at e = 0. 

Average time interval between instantons in the instanton gas for perturbed system can 
be estimated in the following way 

r] « — . (20) 



0J, 



av 



Here u av = ui([Hb + H s ]/2). Near the separatrix following approximation for the frequency 
can be used [0 

""^^^S \H-H,\«H,. (21) 
v ^ m H _ Hs 

Thus interval between instanton transitions is the next 

7TV 1 , \ 

and the density of instantons in instanton gas for the perturbed system can be estimated as 
follows 

r] \2u$ u eis 2 J 

We see that, in contrary to the case of non-perturbed system, for perturbed system the 
density of instanton gas is large if we take into account only exact solutions of Euclidian 
equations of motion. Consideration of approximate solutions only increases the density of 
instanton gas and does not lead to qualitative changes. Comparing of densities for perturbed 
p (p3|) and non-perturbed p (HQ) systems gives 



p _ 1 u 0y /uT o e^ (24) 



Po ~ v / 2 7 n/ + 2u;oln^f 



This ratio is large at large enough ujq (or large enough one-instanton action) and small but 
nonzero e. Thus we obtain that small perturbation (^) can strongly increase the density 
of instanton gas. For calculation of p only exact classical solutions of Euclidian equations 
of motion were taken into account. While to get non-zero value of po we na d to consider 
approximate solutions. Thus the limit e — > can not be applied directly in fl2"4|). Nevertheless 
we have correspondence between results obtained for perturbed and non-perturbed systems 
in this limit. Namely, p (see (fffi)) at e = and po are equal to zero if we take into account 
only exact solutions of Euclidian equations of motion for both cases. 

Let us apply now our formal computation to the physical phenomena in real time. Clas- 
sical instanton solutions in imaginary time (0-0) describe quantum tunneling transitions in 
real time. Some observables can be directly expressed through the instanton density p [|J. 
In particular, rate of the tunnelling between neighbour potential wells (the number of tun- 
nelling transitions per unit of time) and probability of tunnelling are proportional to squared 
density of instantons p 2 . Spectrum and width of the lowest energy zone AE read 

E e w ic; - 2pcos#, < 9 < Ti, AE^4p. (25) 
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Thus squeezing of the instanton gas and increase of the density p in imaginary time mean 
that small perturbation (^) leading to chaos can essentially enhance the tunnelling rate 
and lead to the widening of the energy zone in comparison with non-perturbed system 
(H). Both these results are consequences of the perturbation leading to destruction of the 
single non-perturbed instanton solution and appearance of manyfold of chaotic perturbed 
instantons. An increase of the number of instanton solutions provides a larger number of 
variants for particle to reach one vacuum from another that results in the increase of the 
rate of tunnelling. The decrease of the life-time of the particle in the certain vacuum of the 
system means the widening of the energy zone that is obtained in (p5|). 

On the other hand our results can be considered as a demonstration on the simple model 
of application of the instanton method to the problem of chaos assisted tunnelling. 



4 Conclusion 

In this work we have demonstrated on the example of one- dimensional periodic potential 
that small perturbation leading to chaotic behavior of the system strongly influences on the 
properties of instanton gas. Our estimations show that classical chaos can greatly increase 
the density of instanton gas and rate of instanton tunnelling. 

Both instanton solutions and chaotic behavior can exist in complex systems like field 
theories. The relation between chaos and instantons in such theories is not trivial. Theory 
of strong interactions (QCD) is the most interesting example, where the investigation of 
instanton gas (or instanton liquid) could shed light on the structure of hadrons M, |12[ . 
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